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We compare two advanced designs for gravitational- wave antennas in terms of their ability to de- 
tect two possible gravitational wave sources. Spherical, resonant mass antennas and interferometers 
incorporating resonant sideband extraction (RSE) were modeled using experimentally measurable 
parameters. The signal-to-noise ratio of each detector for a binary neutron star system and a rapidly 
rotating stellar core were calculated. For a range of plausible parameters we found that the advanced 
LIGO interferometer incorporating RSE gave higher signal-to-noise ratios than a spherical detec- 
tor resonant at the same frequency for both sources. Spheres were found to be sensitive to these 
sources at distances beyond our galaxy. Interferometers were sensitive to these sources at far enough 
distances that several events per year would be expected. 



I. INTRODUCTION 



The experimental effort to detect gravitational radiation has advanced substantially since its beginnings in the early 
1960's Two different techniques, resonant-mass antennas and interferometers, have been developed over the years, 
but in the last decade construction has begun on long baseline interferometers designed for very high sensitivity. 

In previous work the sensitivity of these two technologies were compared using models of potential sources of 
gravitational waves. In that work, different diameter spherical, resonant-mass antennas were compared with the 
expected sensitivity of initial LIGO |Q| because it was plausible that both detectors could be in operation in the 
early years of the 2000 decade. Generally, LIGO was found to be more sensitive to these sources, especially the 
inspiraling binary neutron stars. However, at higher frequencies, spheres were shown to provide extra sensitivity 
within a restricted bandwidth. These higher frequencies, about 700 Hz to 5000 Hz, are where the gravitational waves 
from binary neutron star coalescence, rapidly rotating stellar cores, and other sources |13| are found. 

Although gravitational wave sources occur at many different frequencies and amplitudes, binary neutron star inspi- 
rals emit gravitational waves at frequencies accessible to Earth-based interferometers, a few to a few hundred Hertz. 
Neutron star binaries also have an amplitude and event rate that makes it plausible that advanced interferometers will 
detect one or more a year At these frequencies, the neutron stars are many times their own radius apart from 
each other and act as point masses. This makes analytically predicting the waveform possible, but also means that the 
details of the neutron star composition (equation of state, radius, magnetic field, etc.) will not effect the waveform. 
When compact bodies such as neutron stars are close enough together that the gravitational radiation being given off 
carries information about their structure, the frequency is higher (typically above 700 Hz) and the amplitude is lower. 
This is true of the internal motion of compact bodies as well, such as the core collapse of a supernova. Determining 
astrophysically interesting parameters from these higher frequency waveforms will be a major goal of gravitational 
wave astronomy once the first detections have occurred and may require detectors specialized for higher frequency 
response. 

Progress with interferometers and delays with spherical antennas make it more relevant now to compare spheres 
and a more advanced interferometer. One possible upgrade of LIGO would include an additional mirror at the output 
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port, which allows for signal recycling Q or resonant-sideband extraction (RSE) g]. These techniques allow the 
frequency of peak response to be selected over a wide range and the bandwidth of the response to be controlled. 
Such an upgrade would allow the advanced interferometer to operate with a similar strain spectrum as a spherical, 
resonant-mass antenna and have more sensitivity than other interferometer configurations at higher frequencies. 

Here we compare the sensitivity of advanced LIGO, including RSE, with that of eight possible spherical antennas 
with different diameters to determine the effectiveness of each technology as a detector of high frequency gravitational 
radiation. The frequency of peak sensitivity and the bandwidth for the interferometer was chosen to correspond as 
closely as possible to the lowest quadrupole resonance of each of the spheres. Different peak-sensitivity frequencies for 
the interferometer are obtained by varying the position of the signal-extraction mirror (by a fraction of the wavelength 
of the laser light). The fractional bandwidth of a sphere is determined by the choice of transducer, and is typically 
below 20%. To reproduce this narrow bandwidth, the transmittance of the interferometer's signal-extraction mirror 
has been chosen to be relatively low. Once the interferometer has been matched to a given sphere, the sensitivity of 
the sphere and interferometer to different, high frequency, sources can be compared. 

The signal-to-noise ratio for both the spheres and the interferometer was computed using the numerically simulated 
relativistic waveforms of two different gravitational wave (burst) sources: (1) the inspiral and eventual coalescence of 
a binary neutron star system, and (2) a rapidly rotating stellar core undergoing a dynamical instability Both 
of these sources are predicted to have high enough event rates that they could be detected by interferometers within 
the next ten years ||,|-|ll],0 

There are many other sources of gravitational waves that could provide interesting physical and/or astrophysical 
information. The stochastic background of gravitational radiation depends on conditions at the earliest times in the 
universe jl^ and their detection would shed new light on cosmology. Searches for scalar radiation would allow for 
tests of gravity beyond the prediction of general relativity. Both interferometers and spheres may play a role in these 
experiments, but we did not consider these sources in this work. 



II. METHOD 

The spherical antennas were modeled using the same method as described in the previous paper The signal- 
to- noise ratio density was calculated for each sphere using the method of Price jl^ as extended by Stevenson . 
This involved calculating the signal-to-noise ratio density from the strain spectrum of each sphere. The parameters 
that entered the strain spectrum were chosen based, as much as possible, on optimistic extrapolations from values 
demonstrated in operating detectors. 

Eight aluminum spheres were modeled with diameters of 3.25 m, 2.75 m, 2.35 m, 2.00 m, 1.70 m, 1.45 m, 1.25 m, 
and 1.05 m. The 3.25 m sphere weighs 50 tonnes and is the largest solid sphere which can reasonably be manufactured 
and transported. It may be possible to get larger radii spheres, and hence lower resonance frequencies, by building 
hollow spheres [|l^ but we did not consider these detectors. A lower frequency, hollow sphere will have advantages 
with lower frequency sources, especially the inspiral phase of the binary neutron star signal, but not with the higher 
frequency sources we considered here. Each sphere was modeled as having six, three-mode inductive transducers 
arranged in the TIGA geometry |3| which allows for omnidirectional sensitivity. The masses of the intermediate mass 
and the final transducer mass (see Fig. |^) were chosen to give a fractional bandwidth as large a possible for each 
sphere. The transducer was modeled as having a dual, superconducting quantum interference device (SQUID) as the 
first stage amplifier. SQUID amplifiers are currently in use on bar detectors [ p^[pj| ]. A three- mode transducer is being 
developed for use on the Allegro antenna |Q and one has been successfully demonstrated on a test antenna |pl[| . 

The limiting noise in a spherical resonant mass detector comes from two sources; amplifier noise and thermal noise. 
The amplifier noise, which comes primarily from the sensing SQUID, can be separated into additive velocity noise 
and force noise. The additive velocity noise can be written as |^ 

5„(/) = A^„?i27r/o/r„, (1) 

where Nn is the noise number of the SQUID, h is Plank's constant, fo is the resonant frequency of the sphere, and 
r„ is the noise resistance of the transducer. We assumed that the SQUID had quantum-limited noise, and hence a 
noise number equal to one. A quantum-limited SQUID suitable for use in a gravitational wave transducer has not yet 
been demonstrated, although the quantum limit has been reached in SQUIDs with low input impedance 122]. A noise 



number of 24 has been reached at 100 mK in a suitable SQUID when cooled on its own |23|. This noise, however, 
was found to increase significantly when placed in a transducer [|4| . 
The noise resistance was calculated from 

r„ = kAcou/ (47r/o) , (2) 
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where n is an experimentally determined spring constant density given in Table ^ and ^coii is the area 

of the pick-up coils. The noise resistance limits the bandwidth of a sphere with three-mode transducers (when the 
masses are properly chosen) according to jlSt 

(5bw ~ {rnl {2T:foms))^^^ , (3) 

where rUg is the effective mass of the sphere for impedance calculations. This effective mass is given by p^ , p5| 

rus = 5/6 X rup (4) 

where x — 0.301 |^^, mp is the physical mass of the sphere, and the factor 5/6 is appropriate for six transducers in 
the TIGA arrangement p^ . The amplifier velocity noise is shown in Fig. ^ graphed as a strain spectrum. 

The force noise from the SQUID is ultimately detected as motion, and therefore must be converted by the mechanical 
transfer function of the antenna. The output noise can be written 

5/,out(/) =iV„n2^/or„|y22(/)P, (5) 
where the admittance matrix element ^22 (/) for a sphere with three-mode transducers can be written 

2/22 (/) = -27r/« (ci +C2 + Cs~ 47r^/^ (c2 (ci + c^) mi + (ci 4- C2) CsTUs) + 

IGir^ f'^CiC2Csmims) / (—1 + 47r^/^ (cimi + c^toi + Cim2 + C2TO2+ 

Csm2 + CslTls) + leTr"*/** (C1C2TO1TO2 - C2CsTOlTO2 " CiCsmimg- 

CiCsm2ms - C2Csm2ms) + 647r^/^ciC2CsTOirn2TOs) ■ (6) 

In the above, ci is the reciprocal spring constant ki for the spring separating the sphere from the intermediate mass 
of the transducer, C2 is the reciprocal spring constant k2 for the spring separating the intermediate mass from the 
transducer mass, Cg is the reciprocal spring constant kg of the effective spring separating the effective mass of the 
sphere from mechanical ground, toi is the mass of the intermediate mass of the transducer, 7712 is the mass of the 
transducer mass, and rris is the effective mass of the sphere for the lowest quadrupole mode. The masses and springs 
in the transducer are shown in Fig. |l|. The spring constants can be written for each stage in the transducer (j=l,2, 
or s) 

k, - (2^/0)' m, + ^ (27r/) (27r/o) m,/Q,, (7) 

where Qj is the quality factor of the appropriate stage of the transducer and sphere. The Q's depend on composition, 
temperature, and the connections between the masses and springs. The amplifier force noise is shown in Fig. ^ graphed 
as a strain spectrum . 

The thermal noise of the sphere can be written [^j 

S-sph^therm = '^keT^ (j/22 (/)) , (8) 

where T is the physical temperature of the sphere. We modeled the antenna as having a T of 50 mK, although 95 mK 
is the lowest a bar has been cooled at equilibrium |^,^. The term 5R (1/22 (/)) depends on the Q's, with higher Q's 
resulting in lower thermal noise. The sphere and intermediate mass were modeled in aluminum and the transducer 
mass was modeled in niobium. Each of the mechanical Q's was modeled as 40 x 10® [p8|-p0[. Depending on the 
design of the transducer, the final Q can be degraded by the addition of loss from electrical coupling to the SQUID 
circuit [ pT[ |. The sphere's thermal noise is shown in Fig. |^ graphed as a strain spectrum. 
All of the noise sources for the sphere were then combined to determine the total noise; 

'S'tot (/) — Su (/) + 5*/, out (/) + 5'sph, therm- (9) 

This total noise is shown, along with each component, in Fig. ^ graphed as a strain spectrum. 

The gravitational wave signal is applied as force on the spherical antenna, but is read out as velocity of the transducer 
mass. Thus, before a signal can be compared to these noise sources, the gravitational strain must be converted to a 
force and be passed through the admittance matrix of the sphere-transducer system. The comparable signal can be 
written 

5^(/) = ^^^^/'l2/2i(/)M/)P, (10) 
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where Y is Young's modulus for the sphere material, p is the density of the sphere material, n is the reduced cross 
section of the sphere and equals 0.215 for the lowest quadrupole mode ^5|] , h (/) is the frequency-domain amplitude 
of the gravitational wave, and the admittance matrix element 1/21 (/) can be written 

2/21 (/) = 2nfcsi I (-1 + 47r^/^ (cimi + c^mi + Cim2 + C2m2 + 0^777,2 + Csr77s) 

— IGtT*/'* (ciC2777i7772 + C2Cs777 l7772 + ClCsr77l 777s + ClCsr772r77s + C2Cs 7772777s) 

+647r^/^ClC2Cs777 l7772r77s) . (11) 

Table ^ shows all the parameters used in the sphere model. 

The interferometer was modeled using a slightly modified version of the BENCH program |3|]. The noise in the 
interferometer is dominated by three types of noise; seismic, thermal, and optical readout noise. Each noise source 
was modeled using parameters from the advanced LIGO white paper js^], which is scheduled to be implemented in 
2005. The corresponding advanced interferometer is scheduled to begin taking data in 2007. A schematic drawing of 
LIGO with RSE is shown in Fig. | 

Seismic noise is expected to dominate the advanced LIGO noise budget at low frequencies. To reduce the effect 
of seismic noise, each element of the interferometer will be supported by a four-stage suspension which in turn is 
supported from a vibration isolation stack. This vibration isolation will consist of two stages of six-degree of freedom 
isolation. It will use a combination of active and passive isolation with an external hydraulic actuation stage. The 
isolation was designed to make seismic noise negligible compared to other noise sources above some /seismic, expected 
to be 10 Hz. In the model, seismic noise was made extremely high below 10 Hz and vanishingly small above this 
frequency. 

Thermal noise will be the dominant noise source in advanced LIGO in the intermediate frequency band above 10 Hz. 
This noise can be divided into two types; thermal noise from the internal degrees of freedom of the interferometer 
mirrors, and thermal noise from the suspension that supports the mirrors. The mirrors are planned to be made of 
m-axis sapphire, 28 cm in diameter and 30 kg in mass. Sapphire has been found to have much lower internal friction 
than fused silica |3^.35], which is used in initial LIGO. However, sapphire suffers from much higher thermoelastic 
damping than silica |36|. 

The internal mode thermal noise from the sapphire mirror comes from structural damping and thermoelastic 
damping. The noise from structural damping can be found from the loss angle by 

5s..(/) = ^^(C,+C2), (12) 

where fcs is Boltzmann's constant, T is the temperature, / is the frequency, and L is the interferometer arm length. 
The constants C\ and C2 are the overlap between the normal modes of the mirrors and the gaussian-profile laser, 
which has a width wx at the input mirror and W2 at the end mirror. They are found from p7^,B8[ 



^ _ (1 - / °^P (~ (C^^./(2^))') (1 - exp (-4C./7/r) + 4C./7/r exp (-2Gfe/r)) . 

' TrrF ^1 C.Jo(G)' (1 - exp (~2C./7/r))'- 4 (C./iA)' exp (-20/7/0 ' 



^2 / JX,exp(-(C;u;,/(2r))V2) 
' /7Vr^ + 127r/7V^. ^ ^ 



72(1-^)^—^ V 



CJo(C.) 

where a is the Poisson ratio of the mirror material, Y is it Young's modulus, r is the radius of the mirror, I is the 
thickness of the mirror, w is the Gaussian beam width of the laser at the mirror, Q are the zeros of the first order 
Bessel function Ji, and Jq is the zeroth order Bessel function. The value of used is the lowest value measured for 
a piece of sapphire [Q . The thermal noise effects of making a sapphire piece into a mirror are under study, but the 
polishing and especially coating of the mirror are expected to cause some excess loss ps] , ^ . 

Thermoelastic damping also contributes to thermal noise from the mirrors. It is found, in the limit of large mirror 
diameter, from 

Sth if) = . ' ^ [l wl + 1 wl) , 14 

\tt f LCv p J V 
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where a is the thermal expansion coefficient, Cy is the heat capacity at constant volume, p is the density, and k is 
the thermal conductivity. Fused silica is available as a back up material which does not have as much thermoelastic 
loss and has recently been shown to have a (/) as low as 1.8 x 10^® in certain circumstances p| . 

Thermal noise from the suspension, which supports the mirrors below the vibration isolation stack, will be reduced 
in advanced LIGO by replacing the steel slings with fused silica ribbons. Fused silica has much less internal fr-iction 
than steel although with ribbon geometry surface loss limits the achievable dissipation Thermal 

noise from a ribbon suspension with surface loss has recently been considered p5| and the results give thermal noise, 
expressed as gravitational wave stress squared per Hertz, as 



(/) = 32fcBT0dii.9 / (^L,u,m 27r/ (^((271/)' - c^^^J + co^jlnjj , (15) 

where g is the acceleration due to gravity, Lsus is the length of the suspension, m is the mass of the mirror, Wpcn is 
the angular frequency of the pendulum mode, and c/jdii is the diluted loss angle. This diluted loss angle is defined as 



0dii = VY /(127f^is\s) d (0th + 0int) , (16) 

where Y is Young's modulus for the ribbon material, 7 is the ratio of stress in the ribbon to its breaking stress, a is 
Poisson's ratio for the ribbon material, d is the ribbon thickness, 0th is the loss angle due to thermoelastic damping, 
and 0int is the loss angle due to internal friction in the ribbon. Thermoelastic damping in ribbons is found from |46| ] 

0th = 1 , n2 2 ' 

O l + (27r/) 

where a is the thermal expansion coefficient of the ribbon material, C is the heat capacity per unit volume, and Td is 
the time constant for thermal diffusion which in ribbons is given by 

Td = dV (^'^) , (18) 
with D being the thermal diffusion coefficient for the ribbon material. The internal friction in a thin ribbon is given 

by [H 

0int = 0bulk (1 + 64/d), (19) 

where 0buik is the loss angle in the bulk of the ribbon material, and ds is the dissipation depth that characterizes 
the excess loss arising from the surface of the ribbon. The numbers for 0buik and ds in Table H represent possibly 
achievable values, lower values for both have been observed |^,^. Determining realizable values for these parameters 
in advanced LIGO is an area of intense research. 

Optical readout noise in the interferometer can be evident at any frequency in the LIGO detection band. This noise 
source has two separate components: radiation pressure noise from the pressure exerted on the mirrors by the laser 
and shot noise from the inherent granularity (photons) of the laser light. These two noise sources are complementary 
to each other, both depend on the laser power. Recently, optical readout noise in a signal-recycled interferometer has 
been considered from a fully quantum mechanical perspective |4^ . The noise spectrum does differ from the one we 
calculate here, but the difference at high frequencies in a narrowband configuration are negligible. 

The optical power stored in the interferometer is an important parameter for the optical readout noise. There are 
a number of optical cavities in LIGO formed by the different mirrors (input mirrors, end mirrors, power recycling 
mirror, signal recycling mirror, etc.) and each one stores a different amount of power. It is convenient to quote a 
single power, the power incident on the beam splitter, and then calculate the power in different cavities in terms of 
this single value. The power at the beam splitter is proportional to the power out of the laser, P, through the power 
recycling factor 

Pbs = GprP, (20) 

which is found from 

Gpr = 1/ (2iV/3 + ttBs) , (21) 

where aes is the fractional power loss at the beam splitter, f3 is the fractional power loss at each mirror, and A'' is the 
number of bounces that the light makes in each arm, on average. In Fabry-Perot cavities, in the large finesse limit, 
the value can be found from the finesse. 
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N = 2TI-K, 



(22) 



where the finesse T is found from the amphtude transmittance of the input mirror, t\. 

T = 2Ti/{tl + 2l3). (23) 

These equations taken together with the parameters in Table ^ give the power at the beam splitter, 

Pbs = 9.3 kW. (24) 

This power must be kept from being too high because absorption of light in the transmitting mirrors, beam splitter, 
and coatings can lead to thermal lensing. The acceptable thermal lensing limit can be calculated from (including a 
factor of 2 safety margin) |49j 

p ~ ^ ^ ('251 

~ 2dn/dT 1.43AgtBS + l-iAg£ + \Na,o.t ' 

where k is the thermal conductivity of the substrate material, dn/dT is the change in index of refraction of the 
substrate with temperature, Ag is the optical absorption of the substrate, and Ocoat is the relative absorption of the 
optical coating. Using the numbers in Table ||, this maximum allowable power is 

Pmax = 750 W. (26) 

In order to realize the higher power in Eq. (p^), a correction scheme must be utilized that increases Pmax by a factor 
of at least 12.4 for sapphire optics. Research is underway to have such a correction scheme available for advanced 
LIGO ||. 

Both parts of the optical readout noise depend on the response of the coupled cavity system in the interferometer. 
This response can be described by the transfer function between the amplitude of the light in the arm cavity and the 
amplitude of light that enters through the input mirror in each sideband p8| : 

Go.i = 1 - rir2 exp {iTa2nf ) - rir^ exp {i {Ts2nf + S)) + r2rz exp {i ((Ta + Ts) 2nf + S)) , (27) 
Go, 2 = 1 - rir2 exp {-iTa2TTf) - rir^ exp {-i (Ts2iTf + 5)) + r2r3 exp {-i {{Ta + r^) 2nf + S)) (28) 

where ri,r2, and are the amplitude reflection coefficients at the input mirror, the end mirror, and the signal recycling 
mirror, respectively, Ta{= 2L/c) is the light transit time between the input mirror and the end mirror, Ta{= 2Liec/c) 
is the light transit time between the input mirror and the signal recycling mirror with L^ec the length of this signal 
recycling cavity, and 5 is the phase accumulated by the reflected light coming off the signal recycling mirror due to 
its position. The amplitude reflection coefficient for the input mirror can be found from 

rl^l-tl-13. (29) 

The amplitude reffection coefficient for the end mirror can be found from 



= \- (3. (30) 



The amplitude reflection coefficient for the signal recycling mirror, r3, is a tunable parameter as is the accumulated 
phase, 5. 

Radiation pressure noise is largest at low frequencies and, for initial LIGO, is masked by other low frequency noise 
(seismic and suspension thermal noise) ^|. In advanced LIGO, the suspension thermal noise may be low enough that 
radiation pressure is important, but it wi 1 still not be the dominant noise source. Radiation pressure was modeled by 

5,ad (/) = 32Pfc, 27r/A h 4 tl rl (l/|Go,i| + l/|Go,2|)' / ((1 - rirz) {2TrffcmLf , (31) 

where f\ is the frequency of laser light and ts is the amplitude transmittance of the signal recycling mirror. 

Shot noise from the laser is the dominant noise source at high frequencies for intial LIGO, and this will continue 
for advanced LIGO. This noise source was modeled by 

o (r.f f{l-nr2) y^Trfn 

^.hot UJ {f,sm{7Tfra)tlr2h{l/\Go,i\ + l/\Go.,2\)J vPbs ' ^ ' 
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where ry is the quantum efficiency of the photodiode. 

AU of these noise sources were combined to create the total noise curve for advanced LIGO; 



Stot{h) — S'scis + Sint + S'susp + 'S'shot + 'S'rad- (33) 

To produce a numerical estimate of the noise and then a signal-to-noise ratio for a given source, values must be 
provided for all the parameters that go into Eq. (|3^). We used values from the advanced LIGO white paper |Q as 
much as possible. The values chosen for all relevant parameters are shown in Table |^ A graph of advanced LIGO's 
noise compared with spheres is shown in Fig. |^. 

The masses in LIGO are designed to be as close to being in local free fall in the sensitive direction as possible. 
Therefore, the strain from a passing gravitational wave directly gives the change in position of the mirrors. The 
comparable signal, similar to Eq. (|lO|), for an interferometer reads 

= (34) 

Using Eqs. ^ and ^ for the noise of a sphere and interferometer and Eqs. |l^ and |3^ as the comparable signals, the 
signal-to-noise ratio density for each detector can be found from 

a(/) = I](/)/5tot(/). (35) 
Integrating the signal-to-noise ratio density gives the signal-to-noise ratio; 

S/N^/r a{f)df), (36) 



where the angle brackets, (...), denote averaging over gravitational wave polarization and direction. This results in 
a factor of 1/5 for interferometers and a factor of 1 for spheres, as spheres are always optimally oriented. This 
value, S/N, is the figure of merit for a gravitational wave detector and will be used to compare the effectiveness of 
these two different approaches to the different astronomical sources. 

To compare the sensitivities of the the antennas, it is useful to calculate the strain spectral density h (/), 

hif) = a{f)/\hif)\\ (37) 

This quantity cncorporates both the total noise of the detector and the cross section. It is the strain spectrum that 
is shown in Figures ||,^ and |[ 

III. SOURCES 

One category of sources for gravitational radiation at high frequencies (above 700 Hz) is from internal motion of 
compact bodies such as neutron stars. The coalescence and merger of neutron stars as well as neutron star formation 
in collapsing stellar cores are promising sources of detectable gravitational waves. Detecting and analyzing these waves 
promises to teach us a great deal about the physics of strong gravitational fields and extreme states of matter 
Because of the high rotational velocities and strong gravitational fields present in such compact objects, numerical 
simulations must include the effects of general relativity to model the system dynamics realistically enough for use in 
analysis of the data from antennas. To accomplish this goal, a three-dimensional Smoothed Particle Hydrodynamics 
code has been modified to include the General Relativistic "back reaction" . The gravitational radiation from 
these systems is calculated in the quadrupole approximation. 

The component stars of a widely separated binary neutron star system will spiral together due to orbital energy 
losses by gravitational radiation reaction, and eventually coalesce |10 5^,^. Because neutron stars have intense 



self-gravity, as they inspiral they do not gravitationally deform one another significantly until several orbits before 
final coalescence When the binary separation is comparable to a neutron star radius, tidal distortions become 
significant, hydrodynamical effects become dominant, and coalescence occurs in a few orbits. 

The inspiral phase of the evolution comprises the last several thousand binary orbits and covers the frequency 
range / ~ 10 - 1000 Hz. The final coalescence is believed to emit its gravitational waves in the kilohertz frequency 
band range 800 Hz < / < 2500 Hz |9||5^-|59|] . The observation of the inspiral and coalescence waveforms will reveal 
information about the masses and spin angular momenta of the bodies, the initial orbital elements of the system, the 
neutron star radii and hence the equation of state for nuclear matter Pj57|"|59t. 
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Theoretical estimates of the formation rates for binary neutron star systems — with tight enough orbits to merge due 
to gravitational radiation within a Hubble time — can be obtained from empirical rate estimates based on the observed 
sample The most recent study gives a Galactic lower and upper limit of 2 x 10^^ yr~^ and ~ 6 — 10 x 10^^ yr~^, 
respectively . Alternatively, by modeling the evolution of the Galaxy's binary star population, the best estimates 
for coalescence events have been estimated to be as high as 3 x 10"'* coalescences per year in our galaxy, and several 
per year out to a distance of 60 Mpc We have used 15 Mpc as an optimistic estimate and 200 Mpc as a pessimistic 
estimate of the distance antennas will need to look to get about one event per year for this source. 

For the simulation presented here, equal mass component stars are used. Each star is assumed to have a total 
mass of Mt = 1.4 Mq, and equatorial radius, i?oq — 10 km, where Mq is one solar mass. The star is modeled as a 
differentially rotating axisymmetric fluid which use a polytropic equation of state, 

P = kpi (38) 
= fcpi+i/", (39) 

where fc is a constant that measures the specific entropy of the material and n is the polytropic index. The value 
n = 1 is used to simulate cold nuclear neutron star matter. Each star rotates counterclockwise about the z— axis with 
an equatorial surface speed of approximately 0.083 c 0. 

Since the binary system spirals together due to energy and angular momentum losses via the emission of gravitational 
radiation, as the orbit decays, it circularizes radially. Thus, it is a good approximation to assume a circular orbit 
provided the system is expected to have existed for some time [|l]. The component stars used in the simulation 
are initially placed on a sufficiently wide circular orbit (center of mass distance is 40 km) so that tidal effects are 
negligible [ p2| , |63| . Because of the large initial separation, the stars start out in the point mass regime, and as a result, 
their waveforms can be compared directly to the theoretical point-mass inspiral for two neutron stars | |6^ ]. 

Although spin-orbit misalignment in coalescing compact binaries can change the amplitude and modulation of the 
gravitational radiation waveforms, this effect is believed unimportant in the case of a binary neutron star system [ p^ . 
Thus, in the numerical simulation, the orbital and spin angular momentum vectors were assumed to be aligned. If 
the magnetic axis is not aligned with the rotation axis, the cjccta from the coalescence can be trapped within the 
field ||65| . The evolution of the magnetic field configuration during the final inspiral phase of neutron star binaries 
may have significant effects on the frequency and tidal distortion of the coalescence, and hence on the gravitational 



waveforms pq] . The inclusion of this effect is the subject of current research 67 1. 

A series of snapshots of the inspiral and coalescence of the neutron stars along with a graph of the waveform gener- 
ated is shown in Figure ^ The waveform differs noticeably from ones generated with purely Newtonian gravity [|62[ . 
The gravitational wave peak due to the bar formed during coalescence, seen in Fig. ^A)(d), is at a much higher 
frequency; 3700 Hz with General Relativistic back reaction compared to 2500 Hz in the Newtonian case. This peak 
is also broader and less pronounced in the General Relativistic simulation than in the Newtonian. 

Rotational instability during the gravitational collapse of a massive star's degenerate core has long been considered 
an interesting possible source of gravitational radiation. A typical scenario in which such a mechanism can operate 
is the collapse of a rapidly rotating stellar core that has exhausted its nuclear fuel and is prevented from collapsing 
to neutron star size by centrifugal forces. If a significant amount of angular momentum remains in an initially 
axisymmetric core, collapse may be slowed or temporarily stalled by centrifugal forces associated with rotation. If 
the core's rotation is large enough to strongly flatten the core before, or as it reaches neutron-star density, then an 
instability is likely to break the core's axial symmetry The growth of such global rotational instabilities 

provides a means for transport of angular momentum out of the core into the surrounding envelope by transforming 
the core into a bar-like configuration rotating about the short axis, shedding mass in a spiral pattern, thereby allowing 
collapse to continue to a supernova |^,^,^9|-Q . 

The strength of the gravitational waves from such a supernova depends crucially on degree of non-sphericity during 
the collapse, and somewhat on the speed of collapse. If a substantial fraction of the collapsing object's angular mo- 
mentum goes into generating gravitational rather than hydrodynamical waves then the gravitational waves generated 
may be nearly as strong as those generated from a coalescing binary jT^ ]. 

The event rates of Type II Supernovae are fairly well known from observations |l3|] . In our Galaxy, Type II 
Supernovae occur approximately once every 40 years, and several per year out to the distance of the Virgo Cluster 
of Galaxies (about 10 Mpc). Beyond this point, the rate increases roughly as the cube of the distance, where by 
300 Mpc the supernova rate becomes ~ 10^ per year |13 7^|. Although it is unclear what fraction of collapsing 



cores may undergo an instability, if only ~ 1/1000 or 1/10^ do so, this phenomena could be a significant source of 
detectable gravitational radiation . 

For the simulation presented here, the pre-coUapsed stellar core is modeled as a differentially rotating, compressible, 
axisymmetric fiuid which uses a polytropic equation of state. The stellar core is assumed to have collapsed to the 
point where centrifugal hangup occurs, reaching near neutron star densities (polytropic index n — 0.5) prior to the 
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onset of a global dynamical instability. To maximize relativistic effects, the core is assumed to have a total mass of 
Mt = 1.4M0, and equatorial radius, i?oq = 10 km. However, the collapse phase itself is not simulated. The star 
rotates counterclockwise about the axis at at equatorial surface speed of approximately 0.122c 0,^. The event 
was modeled as occurring at a distance of 1 Mpc as an optimistic estimate of the distance necessary for antennas to 
see roughly one event per year A series of snapshots of the evolution along with a graph of the waveform used is 
shown in Figure 0. 

Strongly magnetized neutron stars are expected to form at the end of Type II Supernova collapse. For sufficiently 
strong fields, misalignment between the rotation and magnetic axes can distort the star by trapping the ejecta within 
the field |65|. This can cause a reduction in the angular momentum in a rapidly rotating core through magnetic 
braking, which can remove several orders of magnitude from the initial angular momentum over long enough time 
scales . The inclusion of magnetohydrodynamical effects into the existing numerical models will have significant 
consequences on the stability and subsequent evolution of the post-collapsed object [p7[. 



IV. CALCULATIONS AND RESULTS 



We used Eq. 36 to find signal-to-noise ratios for both spherical resonant mass detectors and the advanced LIGO 
interferometer with RSE detecting the binary neutron star and the rapidly rotating stellar core sources. For each 
source, we calculated seventeen signal-to-noise ratios; one each for a sphere and interferometer configuration at eight 
different frequencies plus one for advanced LIGO in a broadband mode. Although signal-to-noise ratio is ultimately 
the figure of merit for a gravitational wave antenna, the comparison of strain spectra (as in Fig. ^) gives a full 
understanding of the relative merits of the two detectors. We present the signal-to- noise ratio calculations to show 
how each antenna performs astronomically. Comparisons between the two instruments can be done solely on the basis 
of noise performance and cross section. 

The frequencies were set by the choice of diameter for the spheres. The largest sphere, 3.25 m in diameter, has 
a mass of 50 tonnes. The smallest sphere chosen has a diameter of 1.05 m. Any smaller, and the cross section for 
gravitational wave detection (implicit in Eq. p^ ) becomes to small for any realistic sources. The resonant frequencies 
of these spheres are given by H 



/o = 1.62/(7r4ph)VWp): (40) 

where dsph is the diameter of the sphere, Y is the Young's modulus of the sphere material, and p is the density of the 
sphere material. This choice of diameters, then, allows for sensitivity between 795 Hz and 2461 Hz, see Table IH. 

The most sensitive frequency of the interferometer's noise spectrum was adjusted by changing the position of the 
signal recycling mirror. A change in position less than the wavelength of the laser light results in a change in phase 
5 in Eqs. ( p7| ) and (|28|). This, then, changes the frequency characteristic of the shot noise. 

The appropriate 5 was found from setting the derivative of the shot noise with respect to frequency at the resonance 
frequency of the sphere equal to zero. This insures the minimum of the shot noise, which is the dominant noise source 
at frequencies above 500 Hz, will be at the same frequency as the sphere's most sensitive point. Since the frequency 
dependence of the shot noise is all in the function 



|Go(/)|-l/(l/|Go,i(/)l + l/|Go,2(/)|), (41) 



this condition can be written 



= 9|Go {f)\ld6 

= (r3 (ri (1 + rl) sin (2^/t, + 5) - r2 {rf sin (2^/ (r, - r,) + S) + sin (27r/ (r, + r,) + S)))) 
/(2 (1 + rjri + T\rl + r^r^ - ^r^r^ (l + r^) cos (27r/ra) - 2ri (l + r^) cos (27r/r, + (5) 

-f 2r2r2r3 cos (2^/ (r, - r„) + <5) + Ir^r^ cos (27r/ (r, + r,) + . (42) 

Using this equation, the appropriate phase shifts for the signal recycling cavity were found, and are presented in 



Table III, 



The bandwidth of the sphere is determined by the impedance matching between the sphere and the SQUID amplifier. 
With a three stage transducer, this bandwidth is given by 



-'BW 



v/m2/mi, (43) 
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where Sbw is the fractional bandwidth of the sphere in the lossless limit For the choices of masses in the 

transducer from Table 0, this bandwidth becomes 



5bw = 10%. 



(44) 



In the interferometer, the bandwidth of the high frequency response is determined by the reflectivities of the 
input mirror and the signal recycling mirror, ri and ra respectively. The bandwidth of both the sphere and of the 
interferometer was calculated from 



To adjust the bandwidth of LIGO to better approximate the noise spectrum of a sphere, the values ti and ^3 were then 
chosen to get the minimum bandwidth possible. The values used for all frequencies were h = Vom and ts = VomE. 
The bandwidth turns out to be higher than 10%, the sphere's bandwidth, for all peak frequencies. Decreasing leads 
to losses dominating over the transmittance which limits the peak sensitivity. 

The resulting signal-to-noise ratios for the spheres and interferometer configurations are shown in Fig. |^ for the 
binary neutron star inspiral and coalescence, in Fig. ^ for the inspiral and coalescence phases separately, and in Fig. |l^ 
for the rapidly rotating stellar core undergoing a dynamical instability. We also calculated signal-to-noise ratios for 
the sources interacting with LIGO in a broadband configuration optimized for binary neutron star inpiral. This 
involves changing the input transmittance, ti to 0.005, the signal recycling mirror transmittance, to 0.05, and the 
accumulated phase 6 to 0.09. These signal-to-noise ratios are shown in the same figures with a dotted line. 



Interferometers utilizing resonant sideband extraction can be more sensitive than the most sensitive spheres, both 
on and off resonance. This condition remains true even when the effects of random polarization and direction of the 
gravitational wave are accounted for. Spheres are always optimally oriented towards incoming waves. Figures ^ and |l^ 
indicate that this greater sensitivity translates into significantly higher SNR's for the interferometer over spheres for 
the two sources we considered. 

These two figures show how sensitive each technology is to the two sources. A properly sized sphere can detect 
the inspiral signal of a binary neutron star system at a distance of 15 Mpc. This is far enough to reach the nearer 
sections of the Virgo cluster of galaxies. According to optimistics estimates [P,p^ this may be enough to detect one 
event per year. Advanced LIGO can see binary neutron star events out to 200 Mpc with a single detector, the most 
likely distance necessary to see multiple events per year [0 . Advanced LIGO will also be able to see rapidly rotating 
stellar core events at a far enough distance to detect multiple events a year. Depending on the size of the sphere, 
resonant mass technology may also have enough sensitivity to see one or more rapidly rotating stellar core events a 
year as well. 

Detecting the coalescence phase of the binary neutron star event would provide information about the structure of 
these stars (e.g., the equation of state and the equatorial radius). Advanced LIGO properly tuned to a high frequency, 
narrowband configuration provides the highest signal-to-noise ratio for this source. In this mode. Advanced LIGO has 
enough sensitivity to detect the coalescence waveform at a distance of 100 Mpc. This may be enough to actually see 
such an event during the expected lifetime of Advanced LIGO. Advanced LIGO tuned to 1292 Hz, where the highest 
signal-to-noise ratio is obtained, is mostly sensitive to the early stages of coalescence. Choosing both a sphere radius 
of 70 cm and a phase, S, for the interferometer of 0.04613, allows these antennas to be tuned to the 3700 Hz of the 
rotating bar peak. The SNR for the sphere at this frequency is only 4.5 x 10^^ at 15 Mpc. Advanced LIGO has a 
SNR of 6.39 at this distance, but event rate predictions are pessimistic about a coalescence happening this close. 

A comparison of the signal-to-noise ratios found in the previous paper ||^ for binary neutron star events with 
Figures || and |l^ shows that the addition of the gravitational wave back reaction to the model does change the waveform 
of the coalescence phase of the binary neutron star evolution. It is important for deciding the best configuration of 
advanced LIGO to know the details of the coalescence waveform. Other effects, notably inclusion of the magnetic 
fields in the neutron stars and post-Newtonian corrections may change all these waveforms, especially for the 
coalescence and the rapidly rotating stellar core events. We have used the best available models to predict the 
gravitational radiation but further improvements are probable and our results should be seen in this light. 

Despite the sensitivity advantages of interferometers, spheres do have benefits which should allow them to find a 
niche in the field of gravitational wave detection. Having simultaneous detection of a single event by two completely 
different technologies will help confirm signals with marginal SNR's; a near certainty in the early years of gravitational 




(45) 



V. CONCLUSIONS 
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wave astronomy. Having a sphere near to an interferometer site will also allow for correlated searches for stochastic 
backgrounds of cosmological gravitational waves ||75| , |76|] . It is conceivable that such pairing may occur in Louisiana 
between LIGO and Louisiana State University and in Italy between Virgo and an Italian sphere. 

Spheres may be particularly well suited for detecting scalar radiation |]77| , [79| , ^ because of their symmetry properties. 
This would allow for exploration of gravity beyond the predictions of general relativity. The comparatively low cost 
of spherical antennas in relation to interferometers could allow for construction of more individual detectors which 
are located more widely around the globe. The decades of experience working with bar detectors will provide useful 
background for sphere projects. Operation in conjunction with the interferometer network, an array of spherical 
detectors will be a valuable asset to the worldwide effort to develop gravitational wave astronomy. 
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TABLE I. Parameters used in model of spherical, resonant mass, gravitational wave antennas. 



Parameter 


Name 


Value 


Source 


Qs 


Sphere quality factor, Al 


40 X 10® 








Qi 


Intermediate mass quality factor, Al 


40 X 10® 




i 




Q2 


Transducer mass quality factor, Nb 


40 X 10® 








T 


Temperature 


50 mK 




~ 




r„ 


SQUID noise number 


1 








dc 


Sensing coil diameter 


9 cm 






k 


K 


Electrical spring constant per area 


3.78 X 10* N/m^ 








ma/mi 


Mass ratio 


100 








mi/ 1712 


Mass ratio 


100 
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TABLE II. Parameters used for model of interferomctric gravitational wave detector. 



Parameter 
L 

Lycc 
^seismic 

A 
P 

V 

Wl 
W2 

a 

P 



Lsus 

d 

^bulk 
ds 



Name 

Interferometer arm length 
Recycling cavity length 
Seismic noise cutoff 
Wavelength of laser light 
Laser Power 

Photodiode quantum efficiency 

Gaussian width of laser at input mirror 

Gaussian width of laser at end mirror 

Relative power loss in beamsplitter 

Relative power loss at each mirror 

Relative absorption of coating at input mass 

Beamsplitter material absorption coefficient 

Thickness of beamsplitter 

Power transmittance of input mirror 

Power transmittance of signal recycling mirror 

Radius of mirrors 

Thickness of mirrors 

Mirror material loss angle 

Length of suspension 

Total ribbon thickness 

Loss angle for the bulk ribbon material 

Dissipation depth of ribbon material 



Value 

4000 m 
10 m 
10 Hz 
1.064 /xm 
125 W 
0.9 
6 cm 
6 cm 
3.5 X 10"^ 
3.75 X lO"'^ 
1 X 10"^ 
40 X 10"" m"^ 
12 cm 
0.03 
0.005 
14 cm 
12 cm 
5.0 X 10"^ 
0.588 m 
1.7 mm 
3.3 X 10"** 
182 ^m 



Source 



[3311 
[531 



TABLE III. Parameters of the signal recycling mirror to simulate the frequency response of spheres. The transmittance 
of the input mirrors was held constant at = 0.03. The transmittance of the signal recycling mirror was held constant at 
t\ = 0.005. Note that the resonance frequency of the 1.25 m sphere in was a typographical error, the value listed here is 
correct. 



Diameter dsph (m) 


Frequency /o (Hz) 


Bandwidth A///o 


Phase 5 


3.25 


795 


0.170 


0.2271 


2.75 


940 


0.172 


0.1921 


2.35 


1100 


0.182 


0.1641 


2.00 


1292 


0.200 


0.1395 


1.70 


1520 


0.225 


0.1182 


1.45 


1782 


0.254 


0.1005 


1.25 


2067 


0.290 


0.08619 


1.05 


2461 


0.330 


0.07179 
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FIG. 1. Schematic drawing of a spherical, resonant mass, gravitational wave detector with a three-mode transducer attached. 
The sphere mass, ms, is connected to mechanical ground, here the center of mass of the sphere. The gravitational wave acts 
as a force, F, between the sphere mass and ground. Each of the complex spring constants, ks, fci, and k2, includes dissipation 
which gives rise to thermal noise. The transducer mass is next to a superconducting pick up coil which stores a persistent 
current. This current is shunted to the input coil of the SQUID in proportion to the motion of the mass. The SQUID amplifies 
the signal and also serves as a source of wideband noise. 



Ground 




13 



FIG. 2. Strain spectra for a 3.25 m diameter spherical, resonant mass antenna including the components of the noise. The 
dashed line shows the thermal noise at 50 mK, the dotted line shows the velocity (forward action) amplifier noise from the 
quantum limited SQUID, and the dashed-dottcd lino shows the force (back-action) amplifier noise from the same SQUID. The 
solid line is the total noise from the spherical antenna. 
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\felocitv Noise 
Force Noise 



600 700 800 900 1000 1100 

f (Hz) 



14 



FIG. 3. Strain spectra for narrowband interferometers and spherical, resonant mass antennas. 

(A) Both interferometer and sphere are maximally sensitive at 795 Hz, corresponding to the quadrupole mode of an aluminum 
sphere with diameter 3.25 m and a phase shift 5 — 0.2271 for the interferometer. 

(B) Strain spectra for four narrowband interferometers, sensitive at 795 Hz, 1100 Hz, 1520 Hz, and 2067 Hz. Also shown are the 
strain spectra for the four spherical, resonant mass antennas with the same resonance frequencies. The spheres are less sensitive 
than the interferometers at the resonance point, but have roughly the same sensitivity as the interferometers off-resonance. 

(A) 




600 700 800 900 1000 1100 
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1000 2000 3000 
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FIG. 4. Schematic drawing of an interferometric gravitational wave detector equipped for resonant sideband extraction. The 
laser creates a light beam which is sent through the power recycling mirror into two arms of a Michelson interferometer. The 
signal exits at the output port through the signal recycling mirror, which forms an additional Fabry-Perot cavity with the 
Michelson interferometer. Finally, the signal is detected past the signal recycling mirror with a photodiode. 



I I End Mirror 



Beam 

Splitter 



Laser 



Power 
Recycling 
Mirror L , 



Input Mirror 



Input Mirror 



End Mirror 



Signal 

Recycling 

Mirror 



Photodiode 



Output 



16 



FIG. 5. Strain spectra for an interferometeric gravitational wave detector with resonant sideband extraction showing all the 
components of the noise. The dash-dotted-dotted line shows the shot noise from a 125 W laser, the dashed line shows the 
sapphire mirrors' internal mode thermal noise, and the dashed-dotted line shows the thermal noise from the ribbon suspension, 
and the dotted line shows the radiation pressure. The solid line is the total noise. 

(A) The noise components when the interferometer is in a narrowband configuration tuned to 795 Hz. 

(B) The noise components when the interferometer is in a broadband configuration optimized for binary neutron star inspiral. 
(A) 




f (Hz) 

(B) 
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FIG. 6. Inspiraling and coalescing binary neutron stars. 

(A) Particle positions for each neutron star during the coalescence phase. The stars first fall together, reducing the gravitational 
wave amplitude at twice the instantaneous orbital velocity. Then a bar forms which creates an increased amplitude at twice 
the rotational velocity. 

(B) The frequency domain gravitational waveform averaged over source orientation. Notice the slight dip just above 1000 Hz 
from the in- fall and the peak near 3500 Hz from the bar. Each neutron star was modeled as having a mass of 1.4 Mq, an 
equatorial radius of 10 km, and a distance from the antenna of 15 Mpc. 
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FIG. 7. A rapidly rotating stellar core undergoing a dynamical instability. 

(A) Particle positions for the neutron star during the gravitational wave emission which shows the bar shape that develops 
from the instability. 

(B) The frequency domain gravitational waveform averaged over source orientation. The star was modeled as having a mass of 
1.4 Mq, an equatorial radius of 10 km, and a distance from the antenna of 1 Mpc. 
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FIG. 8. Energy signal-to-noise ratios for binary neutron star inspiral and coalescence. This source was simulated interacting 
with spherical resonant mass antennas (shown with asterisks) and interferometers operating with resonant sideband extraction 
(shown with circles) and a broadband interferometer (shown with a dotted line). 

(A) The binary neutron stars were assumed at a distance of 15 Mpc. This distance is the closest estimated for a single event a 
year. Both detectors have high enough SNR's to reach this distance. 

(B) The binary neutron stars were assumed at a distance of 200 Mpc. This distance is enough for multiple events for year, and 
advanced LIGO with RSE does have a SNR high enough to reach this distance. 
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FIG. 9. Energy signal-to-noise ratios for binary neutron star separated into an inspiral and coalescence phase at 15 Mpc. 
The waveform was divided at the dynamical instabihty frequency, when the neutron stars are about 20 km apart. This source 
was simulated interacting with spherical resonant mass antennas (shown with asterisks) and interferometers operating with 
resonant sideband extraction (shown with circles) and a broadband interferometer (shown with a dotted line). 
(A) The binary neutron star inspiral phase. 
(0^1 The binary neutron star coalescence phase. 
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FIG. 10. Energy signal-to-noise ratios for a rapidly rotating stellar core undergoing a dynamical instability. This source 
was simulated interacting with spherical resonant mass antennas (shown with asterisks) and interferometers operating with 
resonant sideband extraction (shown with circles) . The rapidly rotating core event was assumed at a distance of 1 Mpc. 
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